AC DRIVEN JUMP DISTRIBUTIONS ON A PERIODIC SUBSTRATE 
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A driven Brownian particle (e.g. an adatom on a surface) 
diffusing on a low-viscosity, periodic substrate may execute 
multiple jumps. In the presence of an additional periodic 
drive, the jump lengths and time durations get statistically 
modulated according to a synchronization mechanism remi- 
niscent of asymmetric stochastic resonance. Here, too, bista- 
bility plays a key role, but in a dynamical sense, inasmuch as 
the particle switches between locked and running states. 
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The diffusion of a forced Brownian particle on one- or 
two-dimensional periodic substrates provides an archety- 
pal model of relaxation in condensed phase with appli- 
cation to adsorbates on crystal surfaces [|l|, transport in 
superionic conductors [2], acoustoelectric currents in one- 
dimensional channels [3|, dispersion of particles in opti- 
cal traps H , or rotation of molecules in solids || . The 
diffusive process from this category best studied in the 
current literature H- ]lq] is the Brownian motion in a 
one-dimensional washboard potential 



-jx - lUq sin x + F + £(£) 



(1) 



(in rescaled units), where the force terms on the r.h.s. 
represent, respectively, a viscous damping with constant 
7, a spatially periodic, tilted substrate modeled by the 
potential 



V(x,F)=co 2 (l 



Fx, 



(2) 



and a stationary Gaussian noise with zero mean (£(£)) = 
and autocorrelation function (£(£)£(0)} = 2^kTS(t). 

In the low damping limit, 7 <C ujq, (strictly speaking 
it would suffice that 7 < 1.2u>o o) the Brownian parti- 
cle executes multiple jumps, namely over many a lattice 
spacing a = 2n, both in the forward and in the back- 
ward direction p5|Jlq | ; with increasing F the length- and 
time-distributions of such multiple jumps, denoted, re- 
spectively, by N(X) and N(T X ), change from an almost 



exponential curve (locked regime, F < F% = (A/k^ujo 
M) to a kind of bimodal curve with a slow exponential 
tail (running regime F > F 2 ~ 3. 367^0 PJL4I, see also 
Fig. 0). To a good approximation N(X) and N(T X ) 
are related by the scaling law X = (F/^/)T x Q, where 
F/7 is the average velocity of the particle in the running 
regime ||. 

In a recent paper Talkner et al. jl7| addressed the ef- 
fects of a weak periodic drive on the jump statistics in 
the intermediate regime with kT <C F <C lo 2 (exponen- 
tial hopping regime p3|): An enhancement of the mul- 
tiple jumping probability has been predicted in the adi- 
abatic limit of very low modulation frequencies. Taking 
on their challenge to a full investigation, we have simu- 
lated numerically the Langevin equation (LE) (uv with 
the dependent drive 



F -> F(t) = F + AFcos(nt), 



(3) 



AF > 0. Note that the time origin has been set to zero 
without loss of generality and the modulation amplitude 
AF has been kept small compared to the dc term Fq in 
all simulation runs. 

In the present report we focus on the time distributions 
N(T X ) of the multiple jumps at the crossover between the 
locked and the running regime: A significant synchro- 
nization mechanism takes place reminiscent of stochastic 
resonance (SR) in a bistable system |1§|| . 



At low temperatures kT ^C Wq, and in the exponen- 
tial hopping regime, 7 <C u>q and Fq < F\, jumps are 
thermally activated and, therefore, relatively short; their 
direction (cither forwards or backwards) is controlled by 
the intensity of the dc drive Fq relative to F\ |ll],[l(|] . Un- 
der such conditions, the Brownian particle sits in a locked 
state with low mobility, namely 7/1 = j(x)/Fq ^C 1. On 
raising F above the threshold value F 2 ~ 3. 367^0 (-^2 
depends weakly on the temperature) 7/i jumps abruptly 
up to close to unity, thus signalling that the particle 
has switched into the running state with average velocity 
(x) = Fq/j. As the particle flies over the washboard po- 
tential it executes sporadic, short stops at the bottom of 
some potential well. For numerical purposes the diffusing 



particle is deemed as locked (or retrapped) any time it 
sojourns within two adjacent V(x,F) barriers for a time 
interval longer than (27) _1 - the characteristic relaxation 
time of the energy variable. In Fig. U\ we display the nor- 
malized distributions of the jump time durations, N(T X ), 
for different values of Fo, ranging from the locked, well 
into the running regime. As expected, all of our distri- 
butions exhibit exponential tails with time constant T x , 
termed from now on the average jump duration, that in- 
creases exponentially with Fo [[7 15 1. No dramatic change 
in the N(T X ) profile occurs in the vicinity of the thresh- 
olds F\ and F 2 . 

The effect of the weak periodic modulation (g) on the 
jump statistics is illustrated in Fig. g, where N(T X ) is 
plotted for the same ac drive, but increasing values of Fq 
(across the threshold F2). For Fo < F 2 the distribution 
of the jump durations is well reproduced by an expo- 
nential function; for Fq ~ F 2 the curve bends to form a 
broad hump around T x — Tq/2 (here Tq = 27r/$7 is the 
modulation period); finally, for F > F 2 the distribution 
develops a full-blown peak structure, with maxima cen- 
tered at T n = tiTq, with n = 1,2,... superimposed on 
a slowly decaying exponential tail. On a closer look one 
recognizes that even for Fo — Fq, the jump distribution 
is multi-peaked; however, the first peak dominates ap- 
parently over the others, thus indicating a strong drive- 
jump synchronization J19| . More remarkably, the posi- 
tions of the N(T X ) maxima shift from T n = (n — l/2)Th 
at F - F 2 to nT n for F » F 2 |0j. 

The qualitative interpretation of these results is simple. 
For Fo < F2 the length and, therefore, the time duration 
of the multiple jumps is relatively short, - we recall that 
kT <C u)q, - anyway, shorter that the modulation period 
Tq. We are dealing with an adiabatic situation, where 
the slowly varying external parameter F(t) during the 
jump process can be regarded as constant; the ensuing 
distribution N(T X ; t) is expectedly exponential [[12| at any 
times, so that N(T X ) = (N(T x ,t)), (. . .) denoting the av- 
erage over a modulation cycle. To this purpose, it should 
be remembered that the jump rate, namely the number 
of jumps per unit of time, is maximum at Fq ~ F 2 and 
decays exponentially with Fo above threshold |L5|. For 
Fo ~ F 2 the system operates at the unlocking threshold; 
a small ac drive amplitude AF suffices to modulate the 
system, alternately, from the locked to the running state 
and vice versa. No surprise that most jumps last the or- 
der of one half forcing period, hence the appearance of 
a bump in the T x distribution at Tq/2. For Fo > F 2 
the Brownian particle is in the running state most of the 
time, so that only infrequent, very long, jumps may take 
place. Clearly, chances are that rare relocking events 
may happen mostly as F(t) hits a minimum Fo — AF. 
As a consequence, we expect that T X are kind of "multi- 
ple" of Tq (jump-drive phase- locking) . Notice that at the 



threshold the average residence time in the locked state 
T r equals the running time T x , whereas in the running 
regime T r becomes negligible M. 

The dependence of the jump synchronization mecha- 
nism on the amplitude AF of the ac drive is displayed 
in Fig. a. After subtracting the exponential background 
due to the random jumps, we concluded that in the weak 
drive regime, AF <C F), the fraction of the synchronized 
jumps (of any duration T n ) grows quadratically with AF. 
As a matter of fact, changing the ac drive amplitude af- 
fects the background decay constant, as well. At finite 
AF values the background seems to fall off faster than in 
the absence of modulation. This suggests that increasing 
AF has a twofold effect: It modulates more and more 
effectively the jump lengths and, simultaneously, dimin- 
ishes the effective dc drive (Fo , to be defined in the follow- 
ing) below the input value F) by an amount that depends 
on the forcing period. 

Finally, in Fig. ^ we show the dependence of the jump 
distribution N(T X ) on the modulation period. Note that 
all distributions plotted here have been normalized to 
unity after rescaling T x to T x /Tq. The overall f2 de- 
pendence of the synchronization mechanism is appar- 
ent: For short modulation periods the jump-drive phase- 
locking is hardly detectable; on increasing Tq the N(T X ) 
peak structure reaches first its maximum and, eventually, 
merges into the background for exceedingly long forcing 
periods. The qualitative interpretation of these results is 
also straightforward. The case of large Tq corresponds 
to the adiabatic situation illustrated above, when we dis- 
cussed the distributions of Fig. with Fo < F 2 ; the 
case of small Tq values represents the opposite situation, 
where the external time-dependent parameter F(t) varies 
so fast in time, that the diffusing particle only responds 
to its effective average value Fo insensitive to its tempo- 
ral modulation - technically speaking, also an adiabatic 
limit; in the intermediate regime Tn — T x the synchro- 
nization mechanism becomes more efficient, due to the 
interplay between externally driven and thermally acti- 
vated locked-to-running transitions |18| . The apparent 
increase of the background steepness with Tq is a graphic 
artifact introduced by rescaling the units of T x ; when 
plotted versus molecular dynamics time-units all curves 
of Fig. decay with time constants corresponding to an 
effective tilt F such that F - AF < F < F . 

We attempt now at interpreting the numerical results 
reported so far as a manifestation of stochastic resonance 
(SR), - an intriguing phenomenon of nonlinear physics, 
where an optimal amount of noise has the capability of 
enhancing the rate of synchronous switches between the 
local minima of a (possibly, asymmetric) bistable system 
driven by a weak periodic modulation. Stochastic reso- 
nance rests upon a noise controlled phase-locking mech- 



anism which can be vividly illustrated in term of switch- 
time distributions [ [19| . A convenient choice is provided 
by the distributions N\(T±) of the residence times T± of 
the system in either state (denoted by ±, respectively). 
As the modulation forces the switch events by perturbing 
the symmetry of the system, N\(T±) peak respectively at 
T n = t£ 0) + nT n , with n = 1, 2, ... and TJ 0) + t! 0) = T a . 
For a symmetric system, T± = Tq/2. An optimal input- 
output synchronization sets in when the first peak of 
Ni (T± ) at T± dominates over the remaining peaks and 
the exponential random switch background pQ ]. 

In the system at hand all ingredients of the 
synchronization-based (or bona-fide) SR [pll are easily 
recognizable: (i) bistability occurs at the threshold F 2 , 
where the particle mobility switches between two dynam- 
ical states, namely the locked state with 7/j, = and the 
running state with 7/1 = 1. The residence times in the 
two states are denoted here by T r and T x , respectively; 
(ii) the noise control parameter is the equilibrium tem- 
perature T. Note that SR may occur even at low damp- 
ing fL8fl ; (iii) the external periodic drive, AFcos(ftt), Eq. 
(||), tilts the system towards either dynamically stable 
state, alternately, thus determining the multi-peaked T x 
distributions of Figs. |[]j. 

A few well-known SR features are apparent: (1) The 
N{T X ) distributions develop their first (and strongest) 
peak for T n ~ 2T X (see Figs. | and |); (2) The first 
peak, at T x = T x , shifts from close to Tq/2 in the 
symmetric case, Fq = F 2 , to Tq in the running regime, 
Fq > F2 + AF. This is a typical asymmetry effect [ p0[ , 
as on raising Fq above F2, T r gets smaller than T x ; (3) 
The peak structure grows sharper on increasing AF, ac- 
cording to a characteristic SR quadratic law [M . 

Peculiar to our model are the properties following: (4) 
On increasing Fo at AF, T and ft constant, more and 
more peaks shoot up from the background, thus reduc- 
ing the input-output synchronization (see Fig. g). The 
explanation in terms of the SR phenomenology is simple: 
The jump times T x coincide with the residence times in 
the running state, which for F > F 2 gets more stable 
than the locked state; accordingly, T x grows much longer 
than Tq or, equivalently, relative to the directed running- 
to-locked transition, the noise intensity is too low for SR 
to take place; hence the multi-peaked structures plotted 
in Figs. |- I ||||; (5) Contrary to most SR studies, 
here the background constant depends on the ac drive, 
namely on AF and, though to a smaller extent, on Tq. 
A pictorial explanation of such a property is provided 
by Fig. 11.24 of Ref. fq]. Consider the most symmetric 
case with Fq = F 2 ; adding AF to F 2 makes the running 
state more stable than the locked state and, vice versa, 
subtracting AF from F 2 makes the running state less 



stable; however, the imbalance toward the locked state 
at Fo = F2 — AF is not compensated by the imbalance 
towards the running state at Fo = F2+AF; therefore, an 
ac drive contributes to the effective static tilt Fq, with 
Fo < F2, defined implicitly through the average jump 
duration as T X {F , AF, ft) = f x (F ). 
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FIG. 1. Jump distributions in the absence of modulation, 
AF = 0, for different values of the dc drive Fo below (a) and 
in the vicinity of the threshold Fi (b). Fo is expressed in units 
of F2 with F2 — 0.085. Other parameter values: 7/wo = 0.03, 
kT/uo = 0.3. 




FIG. 3. Dependence of the jump distributions on the mod- 
ulation amplitude AF at constant period, Tq — 8 ■ 10 3 , and 
static tilt, F0/F2 = 1.24. Other parameter values are as in 
Fig. 1. 
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FIG. 2. Jump distributions in the presence of an ac 

drive with constant amplitude, AF/F2 — 0.12, and period, 

Tq = 2-rr/Q — 6.4 ■ 10 4 , and for different values of the static 

tilt Fq. Other parameter values are as in Fig. 1. 




FIG. 4. Dependence of the jump distributions on the mod- 
ulation period To. (in units of 10 3 ) at constant amplitude, 
AF/F 2 = 0.12, and static tilt, F /F 2 = 1.24. Other parame- 
ter values are as in Fig. 1. 



